Evolution of the dynamical pairing across the phase diagram of 
a strongly correlated high-temperature superconductor 
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We study the dynamics of the Cooper pairing across the T — phase diagram of the two- 
dimensional Hubbard Model, relevant for high-temperature superconductors, using a cluster exten- 
sion of dynamical mean field theory. We find that the superconducting pairing function evolves from 
an unconventional form in the over-doped region into a more conventional boson-mediated retarded 
form in the under-doped region of the phase diagram. The boson, however, promotes the rise of a 
pseudo-gap in the electron density of states rather than a superconducting gap as in the standard 
theory of superconductivity. We discuss our results in terms of Mott-related phenomena, and we 
show that they can be observed in tunneling experiments. 
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In order to understand high-temperature (high-Tc) su- 
perconductivity one must understand the nature of the 
pairing interaction forming the Cooper pairs. Pairing 
in standard superconductors has been successfully de- 
scribed by the Migdal-Eliashberg (ME) formalism an 
extension of the famous Bardeen-Cooper-Schrieffer the- 
ory (BCS) of superconductivity^. The essential ingredi- 
ent of the ME description is a retarded boson-mediated 
interaction between electrons [ij]. It has been argued, 
however, that high-Tc superconductivity is the product 
of a strongly correlated mechanism H , which gives rise to 
an instantaneous electron pairing[5], radically different 
from the ME description. Recent experimental!!, 0, HJ 



.escript 

and theoretical^, [lfl, studies on high-Tc supercon- 
ductors have focused on this problem, but a number of 
questions remain unanswered. In this paper we enter the 
debate, examining the pairing interaction as one moves 
from the over-doped (OD) region to the under-doped 
(UD) region of the phase diagram of the two-dimensional 
Hubbard Model (a minimal model containing the physics 
of high-Tc superconductors [4J). We look beyond the BCS 
and ME theories, using the cellular dynamical mean field 
theory (CDMFT [H The CDMFT can address the 
full frequency-dependence of the pairing interaction and 
the short-ranged spatial correlation. 

It is well established fl~i| by both experiment [lit and 
theory [I] (including CDMFT studies 0, Ef|) that high- 
Tc superconductors have many properties that are BCS- 
like in the OD region and anomalous in the UD region. 
Little, however, is known about the dynamics of the 
Cooper pairing, it being extremely difficult to identify 
its contribution to physical observablesjl, @. We show 
that in the OD region the pairing function does not dis- 
play BCS or ME forms. Despite the presence of a boson- 
mediated pairing contribution at low frequencies, other 
features are relevant up to an energy scale W* of the 
order of the bandwidth reduced by the strong interac- 
tion. In the UD region, however, the pairing function ac- 
quires an ME form, even if at higher energies (but within 



W*) a pair-breaking contribution appears. We connect 
these findings with the emergence of a pseudo-gap in the 
electron spectra at optimal doping (competing with the 
superconducting gap) and we interpret this in terms of 
Mott-related phenomena (lil 17 1. Finally we show that 
these properties can be observed in scanning tunneling 
microscopy (STM) experiments. The two-dimensional 
Hubbard Model Hamiltonian is: 



(1) 



where Ci„ destroys an electron with spin a on site i and 
jii a = c\ a Ci a . We only consider the nearest-neighbor 
t and next-nearest-neighbor t' = — 0.3i hoppings; fj, is 
the chemical potential. We set the on-site Coulomb re- 
pulsion U — 12t, greater than the bandwidth 8t. In 
CDMFT, H is mapped onto a more easily handled im- 
purity model (2x2 cluster model in our case) of inter- 
acting electrons, embedded in a non-interacting bath of 
fermions, subjected to a self-consistency condition^ . 
We solve the cluster-impurity problem using the Lanczos 
algorithm [l8j], which approximates the non- interacting 
bath with an 8-level parameterization. This method 
works at zero temperature[l9j and provides direct access 
to the information on the real-frequency axis, unlike, for 
example, the quantum-Monte-Carlo method(QMC) (2fj| . 

In this study we consider a paramagnetic transla- 
tionally invariant superconducting state by constraining 
the CDMFT equations. Even if other broken-symmetry 
phases compete for the ground-state (at low doping, an- 
tiferromagnetism is expected to take over, for exam- 
ple), this is a well-defined mean-field procedure which 
allows us to access and study the physics governing 
the superconductivity. We leave open the question of 
which terms could be added to Ti. to make this su- 
perconducting mean-field solution into a real ground- 
state. The output of our CDMFT calculation is the 
frequency-dependent normal G™ t {uj) =<C Cj ff cj CT ^> and 
anomalous Fij(oS) =-C Qt c ii cluster Green's func- 
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FIG. 1: (color online). The superconducting order parameter 
$op (multiplied by t 2 g = 58.9 shown in the inset) and the 
anomalous self-energy value ReE a (c<j = 0) are displayed as a 
function of doping 5 (a). The ReE a (tj) (b) and ImE"([j) (c) 
are displayed on a wide energy range < uj < U /2, from the 
OD region (5 = 0.12, 0.08) to the UD region (5 = 0.04, 0.02). 
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tions and their associated normal Y,™° r (u>) and anoma- 
lous S°"°(u) cluster self-energies — 1, ...,4 in the 
2x2 cluster [IF ' 
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We indicate the Fourier transform 
of the time-ordered ground-state average with <C ... 3>. 
In this paper we focus on the nearest-neighbor anoma- 
lous component S a (cj) (the only one numerically non- 
zero), which directly expresses the pairing function. We 
think that it captures the essential features of the super- 
conducting pairing, which has a d-wave momentum de- 
pendence ~ (cos fee — cos ky) (dominated by the nearest- 
neighbor spatial component [22|), as has been well estab- 
lished experimentally [Tij] and theoretically 3 23|. 

In Fig. [T]a, the superconducting order parameter <& op = 
fToo Fi2(u)du is displayed as a function of doping 5 = 
1 — {rii). <S> op has the expected dome-like shape, with 
a maximum around 5 op t ~ 0.07; this locates the optimal 
doping in our Lanczos-CDMFT calculation (a similar Li 
is also obtained in the QMC-CDMFT results of Ref.0). 



This value is smaller than the 5. 



opt 



0.15 observed in 



cuprate materials, which require more complete models 
to be described in detail. However, as with the cuprates, 
we shall define the region of S > 5 op t[S < S op t] as the OD 
region [UD region] . If we compare Q op with the anoma- 
lous self-energy value T, a (uj = 0) we can already observe 
the unusual non-BCS behavior of the Cooper pairing in 
the OD region. If BCS theory is a good approximation, 
E a (0) ~ g&m,, where g is the strength of the pairing 
interaction[24J. But E a (0) roughly scales like $ op only in 
the UD region (with g ~ 58.9), with no scaling possible 
in the OD region (as is shown by the ^-dependence of g 
in the inset). 

To clarify the observation above we investigated the 
w-dependence of E a (w) (Fig. [Tj3,c). As with previous 
results (il. ITU. E a (w) shows features which extend to high 
energies (ui ~ 6i). Re£ a assumes negative values; these 
are particularly evident in the UD region (around 0.5t < 
uj < 3.5t), indicating a repulsive pair interaction. In 



FIG. 2: (color online). The ratio R (defined in formula [2J) vs. 
uj (a). The low-energy close-up (ui < W*) of ReE a (aj) and 
ImE"(w) is displayed in the OD region (b), close to optimal 
doping 8 pt (c) and in the UD region (d). 



absolute value S a decreases with doping, although Fig. 
QJi shows that ReS a (w — > 0) increases in the UD region 
(in agreement with the QMC-CDMFT results (TJ). 

We now look in detail at the evolution of E a across 
the phase diagram. We first employ Cauchy-Kramers- 
Kroenig relations and define the ratio [9l I22I] 
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R(lu) measures the contribution to the low-energy su- 
perconducting pairing Re£ a (w = 0) obtained from the 
range of frequency up to to. We identify three distinct 
doping-dependent regimes (Fig. [2^). (i) In the OD re- 
gion, R{uj) monotonically increases to uj ~ W* ~ 2t, and 
then decreases, showing a wide hump which overshoots 
one, centered around uj ~ l.5t. The small hump also vis- 
ible (small arrow) around uj ~ 0.5t indicates the presence 
of a boson- like contribution to E a . (ii) Close to optimal 
doping 6 pt, the range of increasing monotonicity of R(uj) 
is reduced to uj < t. For uj > t, R{uj) stays roughly flat 
around one. (iii) In the UD region, a narrow peak over- 
shooting one appears at low frequency uj < t. This form 
of R(oj) is similar to that expected from a boson- mediated 
pairing mechanism 0, 0]. 

R{uj) shows that the pairing builds up for uj < 
W* and hints at the presence of a low-energy boson. 
Bosonic modes have actually been detected in high-Tc 
materials^, 0, @|- Figs. [2jD,c,d, which show the low en- 
ergy part of S a , therefore merit close examination. In 
ME theory the boson exchange between pairing elec- 
trons resonates as uj — > uj~ , a characteristic frequency. 
This fact is made evident by a hump in ReS a (w). For 
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FIG. 3: (color online). The superconducting gap A sc , the 
pseudo-gap A pg , A to t = y Af c + A^ 9 and the characteristic 
boson-frequency tj are displayed as a function of doping S 
(a). The conductance ratio N^ mp ' (b) is a simple ideal case 
(Eq. [3), where only the E a -dependence is relevant. (c) 
and (d) represent more realistic-material case, obtained 
via a periodization[l6l] in momentum space of the CDMFT 
results. 



lu > lu however, the probability of emission of a real bo- 
son increases; this is marked by ReS a (w) changing sign 
and by ImE a (cj) acquiring a peak. This peak occurs at 
il ~ w o + ReA tot (0), ReA tot (0) being the lowest particle- 
excitation available in the superconductor (the supercon- 
ducting gap in a standard BCS superconductor). In the 
OD region (Fig. (2b), S a has the ME form described 
above at low to, but it also has other features at higher 
frequencies. The pair interaction is repulsive for lu > 0.5i, 
where Re£ a changes sign. The Im£ a however has a long 
positive tail until u ~ W* ~ 2t, which means that the 
pairing mechanism involves states in the renormalized 
band of width W*. The superconducting pairing in the 
OD region cannot therefore be described simply within 
the ME theory, which typically shows only one or few 
characteristic boson-frequencies (where Im£ a shows well 
defined peaks). Around 8 op t (Fig. Efc) the long tail in 
ImS a reduces strongly. In the UD region (Fig. the 
tail disappears, and E ano assumes an ME shape (in line 
with the findings of Fig. 1), marking a change in the 
pairing mechanism. 

We now link the evolution of the pairing function 
with the relevant energy-scales of the system. In BCS 
theory there is only one such scale, the superconduct- 
ing gap A sc ~ Uoe" 1 '^ 9 ' (where N Q is the density 
of states at the Fermi level), proportional to the char- 
acteristic boson-frequency to . In our results, however, 
there are two relevant energy-scales 1(| 17 1. The first de- 



rives from Re£ a and represents the superconducting gap 
A sc ~ Zu Re£ a (0) (where Zk nod is the quasiparticle 
residue (lflllj). In our theory, A sc can be measured in 
the regions of momentum space close to "nodal" points 
k„od ~ (±§!±§) by photo-emission(25i l2rjj]. its behav- 
ior is non-monotonic with doping (Fig. (3^-). The sec- 
ond energy scale A pg arises in the UD region from the 
normal component £ raor (see Ref.jlR IttI. [27jj ) : it is con- 
nected with the pseudo-gap observed in the normal state, 
and it is monotonic with 8. It is convenient to define a 
total gap[16| A t ot = (Ap fl + A 2 SC )^ 1 which is measured 
in photoemission[2|| [26| close to the "antinodal" points 
kanod ~ (0, ±7r) or (±7r, 0). In the OD region A sc = A tot , 
but in the UD region the contribution of A pg at the antin- 
odes is most important [3, To clarify whether the 
boson-feature observed in Fig. [5] can be related to any 
of these gaps, we estimate its characteristic frequency 
uj ~ £1 — Atot{k a nod) (at lu — we locate the first maxi- 
mum of ImE°(w) in Fig. [2b,c,d) and plot it as a function 
of doping 8. We find that lu follows the behavior of the 
pseudo-gap A pg (suggestive of the experimental results in 
Ref. 0) rather than the superconducting gap A sc . Com- 
ing from the OD region, lu q is roughly constant, just as 
Apg = 0, while A sc increases. In the UD region, lu 
monotonically increases as doping is reduced; it follows 
the rise of A pg rather than the fall of A sc , in spite of 
the ME form of E a . This is due to the fact that, in the 
UD region, negative values appear in ImE a (oblique ar- 
row in Fig. [2jl) which give a negative contribution to the 
ReE a (0) (Eq. [J) and break the Cooper pairs. 
In our interpretationflF 



171 |27| . these findings are 



Mott-related phenomena. At optimal doping S op t the 
antinodal regions of momentum space undergo a Mott 
transition (A 



mi 



> 0), while the nodal regions retain a 
metallic character. These ideas 11, [28| are reminiscent 
of the orbitally selective Mott transition in multi-band 
models. Close to the Mott transition, antiferromagnetic 
fluctuations are important d, [2^, [3(|, and can mediate 
a pairing channel {c^Cji) to form Cooper pairsQ. The 
proximity to the "antinodal" Mott transition results in 
a pairing process involving electrons on a wide energy 
scale (with Im£ a positive up to lu ~ W*). The high- 
Tc mechanism therefore originates as one approaches a 
Mott transition|4[. Once this latter has taken place at 
S op t, a normal-component particle-hole (c irT Cj a ) channel 
is opened, destroying the Cooper pairing at higher ener- 
gies ( A sc decreases) and favoring the rise of A pg . This is 
marked by the appearance of non-BCS properties in the 
normal component together with a pairing function of a 
more standard ME form. 

Finally, let us look at how the doping-dependent de- 
scription of the pairing function can be measured in ex- 
periments. In ME superconductivity it was possible to 
relate the behavior of the superconducting pairing to 
STM measurements Q of the ratio Nt, between the su- 
perconducting and normal-state tunneling conductances 



4 



dl(u>) I dV , obtained under the same conditions: 



N t (oj) 



N s T mp \Lu) 



dI{uj)/dV\ 
dI{uj)/dV |, 



N sc {uj) 
N nor (u)) 



(3) 



4iJ 



f [Re£ a H 2 -ImS n (w) 2 ] 



where ± ^ „ , . , , , , , 



(uj) — i 5^ fe ImG(/c, is the supercon- 



[normal-state] local density of states. N^ mp ' 



ducting 

is the ideal w — > oo limit, obtained in the simplest 
ME analysis of a d-wave superconductor 0, Q neglect- 
ing the momentum dependence of the band-structure 
and taking N nor (oj) « N nor (0) (a good approximation 
in BCS superconductors). N^ mp ' has the advantage of 
being an explicit function of the superconducting gap 
A sc (o;) ~ Zk nod Y> a (io). Its behavior is shown in Fig. 
[3j3. The sharp drop below unity (due to the Re£ a — » 
and ImS a acquiring a maximum, see Fig. [5]) marks the 
presence of the boson (at u> = O, vertical-up and hor- 
izontal arrows in Fig. ^jp), both in the OD region (in 
agreement with experiments @) and in the UD region. It 
is not possible to see the long tail of the ImS a in the 
OD region, as the real and imaginary parts of S a can- 
cel out. We can however observe two clear features: (i) 
the sudden increase of f2 in passing from the OD region 
to the UD region, mainly due to the increase of w , as 
portrayed in Fig. |3k, and (ii) the rise above unity at 
higher frequencies (marked by a vertical-down arrow) in 
the UD region; this is due to the pair-breaking effect of 
ImE a assuming negative values (in a way opposite to the 
boson-drop, now ImS a — * 0, and Re£ a is maximal). The 
approximations used to extract the simple behavior of 
N^ mp ' no longer hold in the real-case situation, where 
the momentum dependence of the band cannot be ne- 
glected and, above all, N nor (ui) is strongly cj-dependent. 
We performed a full momentum-energy-dependent calcu- 
lation, however, adopting a mixed-periodization scheme 
(introduced in Ref.[16l Il7|). which allows us to extract a 
good approximation for G(k,u) from the cluster results, 
while separating the normal component N nor from the 
superconducting N sc . The resulting Nji for uj > and 
N^p for uj < are displayed in Figs. [3J:,d. The curves 
are now more irregular and asymmetric with respect to 
u) = 0, but the features (i) and (ii) which we have dis- 
cussed for N^ mp ' are still present, showing that it should 
be possible to observe them in experiments @. 

We can conclude that the experiments and theories 
showing unconventional properties in the UD region and 
more conventional properties in the OD region of the 
phase diagram of high-Tc superconductors do not tell the 
entire story. We find an unconventional form of super- 
conducting pairing in the OD region, which evolves into 
a more standard form in the UD region. This crossover 
takes place close to the maximum of the order parameter, 
where a pseudo-gap appears in the one-electron density 
of states. We argue that these phenomena are a con- 



sequence of a Mott-like transition taking place at the 
antinodal points of momentum space. 
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